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Abstract 



Quark light cone momentum distributions in nuclear matter and 
the structure function of a bound nucleon are investigated in the 
framework of the Nambu-Jona-Lasinio model. This framework de- 
scribes the nucleon as a relativistic quark-diquark state, and the nu- 
clear matter equation of state by using the mean field approximation. 
The scalar and vector mean fields in the nuclear medium couple to the 
quarks in the nucleon and their effect on the spin independent nuclear 
structure function is investigated in detail. Special emphasis is placed 
on the important effect of the vector mean field and on a formulation 
which guarantees the validity of the number and momentum sum rules 
from the outset. 



* Correspondence to: W. Bentz, E-mail: bentz@keyaki.cc.u-tokai.ac.jp 
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1 Introduction 



The question of how the properties of hadrons bound in the nuclear medium 
differ from the properties of free hadrons is a very important and active 
field of recent experimental and theoretical research. Examples where such 
medium modifications are important are the nuclear structure functions for 
muon and neutrino scattering [1, 2], proton knock-out processes [3], photo- 
production of pions from nuclei [4], etc. In this work we will be concerned 
with the medium modifications of the spin independent nuclear structure 
functions measured in deep inelastic scattering (DIS) of leptons, that is, the 
EMC effect [1]. It is known [5, 6] that the Fermi motion and binding effects 
on the level of nucleons alone cannot explain the observed reduction of the 
nuclear structure function in the range of Bjorken x between 0.3 and 0.8. As 
a result, one must face the challenging task of describing the nuclear systems 
in terms of nucleons with internal quark structure which responds to the nu- 
clear environment and then investigating whether the binding effects at the 
quark level can account for the observations. 

At the present time a direct description of many nucleon systems within 
QCD itself is out of reach, and so effective quark theories are the best tools 
available. Since we are interested in the effects of the nuclear medium, these 
models should account not only for the properties of the free nucleon, such as 
the structure functions, but also for the saturation properties of nuclear mat- 
ter (NM). Early work on this problem [7] was based on the quark- meson cou- 
pling model. More recently, using the Nambu-Jona-Lasinio (NJL) model [8] 
as an effective chiral quark theory, it has been shown that [9, 10] the quark- 
diquark description of the single nucleon [11, 12], which is based on the 
relativistic Faddeev approach [13], can be combined successfully with the 
mean field description of the NM equation of state (EOS). The completely 
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covariant description of the nucleon facilitates the investigation of off-shell 
effects, and allows a rigorous incorporation of the Ward identities for baryon 
number and momentum conservation from the outset. This is particularly 
important in the description of the DIS structure functions [11]^. The model 
therefore provides a powerful testing ground for investigating the origin of 
the EMC effect in terms of nuclear binding at the quark level. In order to 
achieve such a description of nuclear structure functions, while retaining the 
simplicity of the NJL model description, we will restrict our consideration in 
the present paper to the case of infinite NM. 

The approach which we will follow here is similar in spirit to the previous 
investigations of the EMC effect [7] using the model of Guichon and collabo- 
rators [14, 15], which is based on the MIT bag model description of the single 
nucleon. However, one of the purposes of our present work is to present a 
more systematic study of the effects of the nuclear mean fields, in particular 
of the vector field, on the quark light cone (LC) momentum distributions in 
NM. It is also important to present a framework which is consistent with the 
baryon number and momentum sum rules. We will, as a first step toward 
this aim, restrict ourselves to a valence quark description of the single nu- 
cleon, and to a mean field description of nuclear matter. Thus, for example, 
the effects of the pion cloud [16, 17, 18, 19], which gives rise to the leading 
non-analytic behaviour of the parton distribution functions as a function of 
quark mass [20], will not be discussed here. It has, however, been included 
in the NJL model description of the free nucleon structure functions [11] and 
will be incorporated in the nuclear calculations in future work. 

The outhne of this paper is as follows. In Sect. 2 we discuss the LC mo- 

^The model also gives a clear description of the spontaneous breaking of chiral symme- 
try in normal matter, and of the spontaneous breaking of color symmetry and restoration 
of chiral symmetry in high density quark matter [10]. 
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mentum distributions of quarks in NM by using the convolution formalism. 
Our main task is to explain the important effect of the mean vector field in 
a manner consistent with the sum rules. In Sect. 3 we introduce the model 
used to describe both the single nucleon and the EOS of NM, while in Sect. 
4 we present the numerical results and discussions concerning the origin of 
the EMC effect in our description. A summary and conclusions are given in 
Sect. 5. 

2 Light cone momentum distributions 

The LC momentum distribution of quarks per nucleon in a nucleus with mass 
number A is defined as 

f,/A{xA) = ^e^^— < A,P\m7^i^i^-)\AP> . (2.1) 

Here is the total 4-momentum of the nucleus and P_ its "minus compo- 
nent ^, and i/j is the quark field. Because we consider only the isoscalar case, 
N = Z, without the Coulomb force, we need only the isospin symmetric 
combination, fg/A = fu/A + fd/A- 

Equation (2.1) corresponds to the quark 2-point function in the nu- 
cleus, where the minus component of the quark momentum is fixed as A;_ = 
P-Xa/A, and the other LC components are integrated out. The variable xa 
is therefore A times the fraction of the minus component of the total mo- 
mentum carried by a quark: xa = Ak^/P^. We will work in the rest frame 
of the nucleus, where P^ = {Ma, 0), and 

P_/A = Mn/V2 (2.2) 

^Our notations for LC variables are = (a° ± a^)/V^, a± = (oq ± as)/^. We will 
frequently call p- {p+) the "minus (plus) component" of the 4-vector p^. 
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where Mtv = Ma/ A is the mass per nucleon. In this frame we have xa = 
V2A;_/Mjv, or 

Mm ,^ 
xa==-x, (2.3) 

where x is the ordinary Bjorken variable for DIS on a free nucleon with mass ^ 
Mm = 940 MeV. 

Since in the nuclear medium there is no Lorentz invariance with respect 
to the motion of a single nucleon, it is natural to use a non-covariant nor- 
malization of state vectors. In this work we use the "non-covariant LC nor- 
malization", where 

< A,P|V^7+-0|A^ 3A (2.4) 

In this normalization, the vertex function of a bound state, like the nucleon, 
is defined by the residue of the i-matrix in the corresponding channel at 
the pole in the plus component of the momentum of the bound state - see 
Subsect. 3.1 for the case of the nucleon. 

In this Section we wish to discuss some features of the distribution (2.1) 
which are largely independent of the details of the model for the nucleon and 
NM. There are two approximations which we need for this purpose. The first, 
which corresponds to the mean field approximation, is that in the presence 
of the nuclear scalar and vector mean fields the energy spectrum of a single 
nucleon with 3-momentum p has the form 



Po 



^ep^^fMl^+^ + 3Vo = Ep + 3Vo. (2.5) 



Here the effective nucleon mass, Mjv = Mn{M), is obtained as a function 
of the effective (constituent) quark mass, M, by solving the bound state 
cciuatiou for the nucleoli using the solution of the in-medium gap equation 
^We will use the subscript to denote masses at zero baryon density. 
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for M . We denote by = (Vq, 0) the constant mean vector field acting on 
a quark in tlie nuclear medium. That is, the quark Hamiltonian in the mean 
field approximation for NM at rest has the form 

Hg = hg + VoQ, (2.6) 

where hg is the quark Hamiltonian without the mean vector field, and 
Q = J d^x-'tp''(x)ip(x) is the quark number operator. We also define the "ki- 
netic momentum" of the nucleon in the medium as 

pi^^p^^-SVr (2.7) 

The second assumption is that the distribution (2.1) can be evaluated by 
using the familiar one-dimensional convolution formula 

fq/A{xA) = j dyAj dz6{xA - Vaz) fq/N^z) fN/AiVA) , (2.8) 

where fg/^iz) and fN/AiVA) are the LC momentum distributions of quarks 
in the nucleon, and of nucleons in the nucleus (per nucleon), respectively: 

dw" 



/(i'li; _ 
^ N,p\mi^i^i^')\N,P >, (2.9) 

fN/AiVA) = ^/^e^^-^---/^< AP|«0)7+V'iv(^-)|AP> ■ 

(2.10) 



Here V'at is the nucleon field and the normalizations are 

< N,p\ii-f+'ip\N,p> = 3, (2.11) 
<A,P\i^j,j+iPN\A,P> = ■^<A,P\ipliPr^\A,P>^A, (2.12) 



*In our quark-diquark model for the nucleon, which will be discussed in Subsect.3.1, 
the interpolating field for the nucleon is local - see Ref. [10] for the explicit expression. 
We also note that the main assumption which leads to the convolution formula (2.8) is 
that the LC momentum distribution of a nucleon with virtuality i"^ has a sharp peak at 
£^ = M'^, so that the quark momentum distribution in the nucleon can be evaluated at 
i'^ = and taken outside the integral over the virtuality. 
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where the first equahty in Eq. (2.12) holds in the rest system of NM. It 
follows from the expressions (2.9) and (2.10) that z is the fraction of the 
minus component of the nucleon momentum (p'') carried by a quark, z — 
k^/p-, and ija is A times the fraction of the minus component of the nuclear 
momentum (P'*) carried by a nucleon, = Ap_/P_ = •\/2p_/Mjv. 

In the following two Subsections we will discuss the distributions fq/N{z) 
and fN/AiVA) separately. In particular, we focus on their dependence on the 
mean vector field and the validity of the baryon number and momentum sum 
rules. In Subsect. 2.3 we will return to the convolution formula (2.8). 

2.1 Quark LC momentum distribution in the nucleon 

If one has a model to describe the nucleon as a bound state of quarks, the 

evaluation of the distribution fq/N{z) of Eq.(2.9) can be reduced to a straight- 
forward Feynman diagram calculation by noting that it can be expressed as 
follows [5, 11] 5; 

U/n{z) = -'j^j['-y) k)) , (2.13) 

where the quark two-point function in the nucleon is given by 

N,p), (2.14) 

and the trace in (2.13) refers to the Dirac and isospin indices. In the quark- 
diquark model, which will be introduced in Sect. 3, such a direct calculation 
is based on the Feynman diagrams of Fig.l - as discussed in Appendix A. 

The difference with respect to the free nucleon case considered in Ref. [11] 
is that all propagators have to be replaced by those including the scalar and 
vector mean fields. The effect of the scalar field can easily be incorporated 

^We remind the reader that for connected LC correlation functions the T-product is 
identical to the usual product [21]. 
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M, 



{N,p\T t/j^M^) 




Figure 1: Graphical representation of the function fq/N{z) of Eq.(2.13) in the 
quark- diquark model. The single (double) line denotes the quark propagator 
(scalar diquark t-matrix), the hatched circle stands for the quark-diquark 
vertex function and the operator insertion, denoted by a cross, stands for 
^'^5 {z — k^/p-). 



into the effective masses, but the effect of the vector field is a bit more subtle. 
In this Subsection we wish to determine this dependence on the vector field 
more generally, without explicit reference to Feynman diagrams. 

For this purpose, we note the following two points. First, the quark field 
satisfies 

il^{w) = e^^''-"'^(0)e-^^''-^ , (2.15) 

where Pj^ — {Hg, Pg) the 4- momentum operator for quarks. The dependence 
of the quark Hamiltonian on the vector field has been given in Eq. (2.6). Since 
[hg, Q] — and Q is time independent, we have 

^^^inJ'[Q,^|;{nJ)]^-iw''^^;{w), 

the solution of which can be written in the form of a local gauge transforma- 
tion ^: 

^jj(w) = e-'^-'^i/joiw). (2.16) 
Here ipo is the quark field in the absence of the vector potential. 



^To check the relation (2.16) in classical field theory, it is sufficient to note that 
ipQ {i^ — M) ipo = ^j} {i^ — Y — M)il}, i.e., the local gauge transformation generates the 
coupling of the quarks to the (constant) vector potential. 
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Second, if the nucleoli state \N,p) is an eigenstate of the quark Hamil- 
tonian, Hg, with eigenvalue p° given by Eq.(2.5), and of the quark number 
operator Q with eigenvalue 3, it is also an eigenstate of hg with eigenvalue p^, 
which is related to p° by Eq.(2.7). We can express this fact in the following 
way: 

|iV,p) = |iV,pN)o, (2.17) 

where the subscript characterizes a nucleon state without the external 
vector field, which we normalize in the same way as Eq. (2.11). If we insert 
the relations (2.16) and (2.17) into Eq.(2.9), we obtain 

fM^) = ^ Uim{^' ^^'-^)^ (2-18) 

whcrc fq/No{z') is the distribution function in the absence of the vector po- 
tential and can simply be evaluated by following the method for the free 
nucleon case [11] with modified masses. Since in our valence quark picture 
(no sea quarks) the function fg/m{z') has support for < z' < 1 [11, 21], 
the distribution fg/N{z) has support in the region Rn{z) defined by 

Rn{z) : — < z < 1 - ^— (support of fq/N{z)) . (2.19) 

Here we make several comments. First, in the quark-di quark model, the 
result (2.18) also follows from a straightforward calculation of the Feynman 
diagrams of Fig.l - as shown in Appendix A. Second, this discussion con- 
cerning the dependence on the vector field can easily be generalized to any 
hadron (H) with quark number H. This leads to 

fM^) = ^ h/m{z' =^z-^), (2.20) 
Ph- Ph- Pn- 

where p'^ is the momentum of the hadron and 

p'^=p^'-HV'' (2.21) 
10 



is its kinetic momentum. The distribution (2.20) has support in the region 
Rh{z) defined by 

Rh{z) : ^ < ^ ^ 1- ^^"^^^" (support of fg/H{z)) • (2.22) 

The quark number, H — 3, corresponds to the case of the nucleon (H=N), 
which was discussed above. Third, the relation (2.18) shows that fq/N{z) ac- 
tually depends not only on but also on or on ua (since ua = Ap^/P^), 
although we do not indicate this fact explicitly in our notation. This depen- 
dence on the nucleon momentum, which arises because in the medium there 
is no Lorentz invariance with respect to the motion of a single nucleon, has 
to be taken into account in the convolution integral (2.8), as will be discussed 
in Sect. 3. Fourth, in Appendix B we show the consistency of (2.22) with the 
spectral representation 

/,/^(^)=p_v^^5(p_(l-z)-pL")) (2.23) 

n 

where ■0(+) = -^+'0 with A+ = 7^7+/ V^. That is, following the discus- 
sions in Ref.[21], there are contributions from connected and semi-connected 
diagrams to the r.h.s. of Eq.(2.23) and the connected piece (diagram 1 of 
Fig. 8 in Appendix B) has support only in the region (2.22), while the semi- 
connected pieces have support in different regions of z. The connected piece 
can be interpreted as the LC momentum distribution of quarks in the hadron 
H in the sense of the parton model. 

If we suppose that the distribution fq/m{z) satisfies the number and 
momentum sum rules, which is the case in the quark-diquark model for the 
nucleon [H — N) [11], Eq.(2.20) can be used to confirm the validity of these 
sum rules for fq/H^z): 

f dzfq/H{z) = [' dz' fq/H0{z') ^ H , (2.24) 
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/ dzzfq/H{z) = dz' Z' fq/Ho{z') + — dz'fq/Ho{z') 

JRh P- Jo P- Jo 

= ^ + ^ = 1. (2.25) 
P- P- 

This derivation clearly shows that the shift of the momentum fraction induced 

by the vector potential (see Eq.(2.20)) is necessary to satisfy the number and 
momentum sum rules at the same time. 

2.2 Nucleon LC momentum distribution in nuclear mat- 
ter 

The function fN/AiVA) of Eq.(2.10) has been discussed in detail in Ref.[6] 
for the case of a mean field description of NM and we briefly summarize the 
derivation here for convenience. First we express (2.10) as follows: 

where we used Eq.(2.2) and introduced the nucleon two-point function in the 
medium 

Gn,M -ijd'w e'P-^ {A, P\T [^^_ JO) i^NA^)] \A P) ■ (2-27) 

The relation to the familiar Feynman propagator in the medium Sn{p) [22] 
is given by ^ 

GM{p)^Vy/2SM, (2-28) 

where V is the volume of the system. The propagator S^ip) consists of a 
"Feyumau ])art'" Syj.-{p) = (j^^ — Mjv + and a density dependent part, 



''To verify Eq.(2.28), we note that the usual (dimensionless) NM ground state \p > at 
the saturation density p satisfies (pIV'atV'wIp) = P- Comparison with the normahzation 
(2.12) gives \A,P) = VF2i/4|p), ^]^^^^ ^^g^^g Eq.(2.28). 
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SNnip), which takes into account the fact that NM consists of nucleons with 
momenta up to the Fermi momentum pp: 

Snd{p)=^-^^^17^5{p^o-E,) Q{pf-\p\). (2.29) 

where the kinetic momentum pn in these expressions was defined in Eq.(2.7). 
In the mean field description of the NM ground state, we have to replace [22] 
Sn — > Snd in loop integrals like in Eq.(2.26), which becomes 

fN/AiVA) -^j-^j[y^- ^) {l^S^nip)) . (2.30) 

Here p is the baryon density and we used the fact that at the saturation 
density the mass per nucleon is equivalent to the Fermi energy ^ that is, 
using Eq.(2.5), 



Mt^ = tF = ^pI + MI + Wo = Ef + 3Vo. (2.31) 

The expression (2.30), which is graphically represented by Fig. 2, can easily 
be evaluated with the result 

" " 2 



(2.32) 



l-^<y^<l + ^. (2.33) 

dp 6f 

It is easy to verify that (2.32) satisfies the relation ^ 

fN/AiVA) = ^Jn/ao {v'a = - ^) ' (2-34) 



^This is the "Hugenholtz-van Hove theorem", which follows simply from the definition 
of the chemical potential {ep = d£/dp, where £ is the energy density of NM) and the 
saturation condition for Mj\f = £/ p: 

d £ cp £ ^ — 
^ — = — -^=0^ep = Mn- 
dpp p p^ 

^The relation (2.34) is formally similar to (2.20), which refers to a hadron in an external 
vector field. However, one should not push this analogy too far, since no stable NM state 
exists without its internal vector field and therefore it does not make sense to define the 
quantity fN/AoiUA) simply by setting Vq = directly in Eq.(2.10). 
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Bp 

P 



X 




Figure 2: Graphical representation of the function fN/AiUA) of Eq. (2.30). The 
closed line denotes the nucleon propagator Snd, and the operator insertion, 



denoted by a cross, stands for 7+ (5 (j)^ 



V2 



where fN/AoiUA) distribution obtained by setting Vq — in Eqs. (2.32) 

and (2.33): 



fN/Aoiy'jd 



.Pf, 



Pf 



- (1 - Va)' 



^- -FT <yA<^ + -Er- 

Jllp Jllp 



(2.35) 
(2.36) 



The number and momentum sum rules for the nucleons arc easily verified 



by noting that any function of the form f{y) 



3 1 
4^ 



(1 — yy , where 



1 — a<y<l + a, satisfies J dy f{y) — J dyy f{y) — 1 for any value of a. 
That is, both f^/A and /jv/ao satisfy the baryon number and momentum sum 
rules 



/ 



j dyAfN/A{yA) = j dyAyAfN/A{yA) = 1, 
dy'AfN/Aoiy'A) = dy'^y'^ fN/Ao{yA) = : 



(2.37) 
(2.38) 



where the ranges of integration are given by Eqs. (2.33) and (2.36). 

Although the form (2.32) formally satisfies the sum rules for any values 
of Pi? and ep, one should keep in mind that it is valid only at the saturation 
point - since we used the relation (2.31) which follows from the saturation 
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condition. Neither the momentum sum rule, nor the relation (2.34), would 
hold if one were to take Mjv to be a fixed parameter. 



2.3 Quark momentum distribution in nuclear matter 

In order to obtain the dependence of the quark LC momentum distributions 
in NM on the vector field, we have to insert (2.18) and (2.34) into the con- 
volution formula (2.8). As we pointed out already, because of the presence 
of the mean vector field, which refers to the NM rest frame, fq/N{z) actu- 
ally depends also on y^. Convolution integrals of the type (2.8) are treated 
generally in Appendix C and the result is as follows: 

fq/A{XA) = ^fq/AQ{x'A = - ^) ■ (2.39) 

Hip Hip Hip 

Here the distribution without the explicit effect of the mean vector field is 
denoted by fq/AQ^x'^) and defined as 

fq/A0{x'A) = I^VaI K^'a - V'a^) fq/miz') fN/AoiVA) ■ (2-40) 

The distribution fq/Aoix'A) includes, besides the Fermi motion of nucleons, 
only the effect of the scalar field, which is incorporated easily into the effective 
masses. Because of the constraint < z' < 1 and the support of given in 
Eq. (2.36), the distribution (2.40) has support for < x'^ < 1 +pp/Ep and 
the distribution (2.39), which includes the effect of the mean vector field, has 
support in the region 

Ra{xa) ■ — <xa< ^^^'P^^ (support of fq/A{xA)) ■ 

ep tp 

(2.41) 

Since we have already shown the sum rules for /^/jv and /at/a, the sum 
rules for fq/A{xA) follow immediately from the convolution formula (2.8). One 
can also make direct use of Eq. (2.39) to check the validity of the sum rules, 
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which shows the importance of the shift arising from the mean vector field 
more exphcitly. Indeed, using the fact that the distribution (2.40) satisfies 
the number and momentum sum rules, we obtain 

/• fl+ao 

J^^dxAfg/A{xA) = I dx'^fg/Aoix'A)^^, (2.42) 

/ dXAXAfq/A{xA) = / dx'Ax'Afq/Ao{x'A)-\ / dXA fq/Ao{x'A) 

JRa 6p jo €p Jo 

= ^ + ^ = 1, (2.43) 

ep ep 

, Pf 
where oq = 

Ep 

3 Model for the nucleon and nuclear matter 

In order to evaluate the quark distributions and structure functions in NM, 
we need a model for the nucleon and the EOS of NM. In this work, we use 
the NJL model as an effective quark theory to describe the nucleon as a 
quark-diquark bound state and NM in the mean field approximation. The 
details are explained in Refs. [9, 10], and here we will only briefly summarize 
those points which will be needed for our calculations. 

The NJL model is characterized by a chirally symmetric 4-fermi interac- 
tion between the quarks. Any such interaction can be Fierz symmetrized and 
decomposed into various qq channels [13]. Writing out explicitly only those 
channels which are relevant for our present discussion, we have 

jC^i^ii^-m)iP + G^ (^{i^^jf - {i^M^pf^ - {i^n^y + . . . (3.1) 

where m is the current quark mass. The quark bilinears j/jiI) and ip^^jl) have 
expectation values in the NM ground state \p), which we will separate as usual 
according to ipip = {p\ipilj\p) + (: ipip :) and ■ipy'^ip = (plipj^iplp) + (: ijj'j'^ip :). 
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The Lagrangian can then be expressed as 

£^?»-M-f,*.-i^^ + ^ + £. (3.2) 

where M = m - 2G^(p|V^V|p) and ^ = 2Gu,{p\^p-f^ij\p), and Ci is the 
normal ordered interaction Lagrangian, including the counter terms linear in 
the normal ordered products. 

In order to fully define the model one has to define a cut-off procedure. 
In all calculations in this paper we will use the proper time regularization 
scheme [23, 9], where one evaluates loop integrals over a product of propaga- 
tors by introducing Feynman parameters, performing the Wick rotation and 
replacing the denominator (A) of the loop integral according to 

dTT"-^e-^^ (n>l), (3.3) 

A- (n-l)!7i/A?,v 
where Air and Auv are the infrared (IR) and ultraviolet (UV) cut-offs, re- 
spectively. Although there are no IR divergences in the theory, the IR cut-off 
plays the important role of eliminating the imaginary part of the loop inte- 
grals. That is, it eliminates the unphysical thresholds for the decay of the 
nucleon and mesons into quarks [23] , thereby taking into consideration a par- 
ticular aspect of confinement physics. It has been shown in Ref.[9] that this 
is crucial to describe saturation of the NM binding energy in the NJL model. 

3.1 Model for the nucleon 

One can use a further Fierz transformation to decompose Cj into a sum of qq 
channel interaction terms [13]. For our purposes we need only the interaction 
in the scalar diquark [J^ = 0^,T = 0, color 3) channel: 

jCi,, = Gs {i^ T^^^ij") (^^ (C-^s) T2/5^^) , (3.4) 
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where 13^ = \JS/2 A'^ {A = 2, 5, 7) arc the color 3 matrices and C = ij2lo- 
The couphng constant Gs will be determined so as to reproduce the free 
nucleon mass. 

The reduced t-matrix in the scalar diquark channel is given by [9] 

"''''^ i + 2G.n.fa) ="°'''°' P-'' 

with the scalar qq bubble graph 

n.(g) = 6t J 0^,trn b,S{k)^,S {-{q - k))] - n,o(gD) ■ (3.6) 

In the second equalities of Eqs. (3.5) and (3.6) we indicated that these 
quantities can be obtained from their expressions without the effect of the 
mean vector field by replacing q ^ Qd, where qj) is the kinetic momentum 
of the diquark defined by H — 2 in (2.21). This is easily seen by using the 
following form of the quark propagator in the presence of the mean vector 
field: 

S{k) = (^Q - M + ze) = 5o(fcQ) , (3.7) 

where the kinetic momentum of the quark is defined by if = 1 in (2.21) and 
So is the usual Feynman propagator without the effect of the external vector 
field. 

The relativistic Faddeev equation in the NJL model has been solved nu- 
merically for the free nucleon in earlier work [13]. For the present purposes, 
however, because we need to solve the equation many times in order to obtain 
self-consistency in-medium, we restrict ourselves to the static approximation 
to the Faddeev equation [24] , where the momentum dependence of the quark 
exchange kernel is neglected. The solution for the quark-diquark t-matrix 
then takes the simple analytic form 

18 



with the quark-diquark bubble graph given by 

niv(p) = - / ^ S{k) r{p -k)^ UmiPN). (3.9) 

In these relations, again, the subscript characterizes the quantities without 
the effect of the mean vector field. The nucleon mass Mjv is defined as the 
pole of (3.8) at = Mjv and the positive energy spectrum has the form 
given in Eq.(2.5). 

The nucleon vertex function in the non-covariant LC normalization is 
defined by the pole behavior of the quark-diquark t-matrix in the variable 
p+: 

r.(rtr.(p) ^p^^j (3.10) 

P+ - ep 

where ip is the "LC energy". From this definition and Eq.(3.8) one obtains 



^n{p) = \ -Zn — ^ Miv(pN) , (3.11) 
V Pn- 

where wjv is a free Dirac spinor for mass Mj^ normalized as unUn — 1 and for 
an on-shell nucleon we have pn- = {Ep — ps) / -\/2. The normalization factor 

Z]\f is given by 

= f f-T^) 1 " ■ (3-12) 

We note that in this normalization the vertex function satisfies the relation 

TM (^5^) r„(P) = 1. (3.13) 

In the numerical calculations of this paper, we will approximate the quantity 
To by 

ro(g) - UG, - (3.14) 



^°Here, and in the following, the derivative w.r.t. ^ oi a Dirac matrix valued function 
a(p^) + ■fi){p'^) is evaluated by using = as 2^a'(p^) + 6(p^) + 2p^6'(p^), where the 
prime denotes a derivative w.r.t. p^. 
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Here the diquark mass Mo is the pole of tq of Eq.(3.5) and the residue is 
given by 



3.2 Evaluation of quark distributions via the moments 

The calculation of the distribution fq/m-i which appears in the convolution 
integral (2.40), in the quark-diquark model follows the lines explained in 
Ref.[ll] for the free nucleon. The difference, however, is that here we are 
using the proper time regularization scheme, which is defined in the Euclidean 
formulation after Wick rotation. In this case, the regularized LC momentum 
distribution cannot be calculated directly from the Feynman diagrams of 
Fig.l, since in the loop integrals one has to fix the minus component of the 
quark momentum to a real value, see Eq.(2.13), which is not possible in the 
Euchdean formulation. We therefore consider first the moments 



and then derive an expression for the quark distributions which can easily 

be used in Euclidean regularization schemes^ ^. 

In the following, we will outhne the calculations for the "quark diagram" 

(first diagram of Fig. 1) and the "diquark diagram" (second diagram of Fig. 

1) separately, leaving the details to Appendix D The contribution from 

the quark (diquark) diagram will be denoted as fg^^oi^) ifq^mi^)) ^^'^ 

^^A more general discussion of this method can be found in Ref.[26]. 
Since all formulae given below in this Subsection refer to the case without the explicit 

effect of the mean vector field, there is no difference between the kinetic and canonical 
momenta and therefore we denote the nucleon, quark, and diquark momenta as simply 
p, k, and q, respectively. 




(3.15) 




(3.16) 
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total isoscalar quark distribution in the nucleon is then given by the sum 

f,/m{x) = + fj^^Ux) . (3.17) 

Here, as in the previous subsection, the subindex (0) indicates that the effect 
of the mean vector field is not included. The contribution of the quark 

diagram to the n-th moment of fq/m{x) is given by 

A^^' = -riv(p)^ (/ (04 (^)" ' '^o(^) ro(p - k)^ T^ip). (3.18) 

From (3.13) and (3.9) it follows that the contribution of the quark diagram 
to the first moment (quark number) is equal to 1. 

Because of the derivative in (3.18), only the pole term in the expres- 
sion (3.14) for To contributes. Introducing a Feynman parameter (a), and 
observing that g = 0, we obtain (see Appendix D) 

J {2t:Y Jo {k^-A{a,p'')f 

(3.19) 

where for the calculation of the derivative it is understood that = and 
at the end one sets^ = M^. The quantity A^ajp"^) is defined as 

A{a,p^) =M^{l-a) + Mla-p'^a{l-a). (3.20) 

From (3.12) we obtain the following form of the normalization factor Zn' 

Z7} = -ioD TT, / 7 — 77 / da '-^ , 3.21 

^ J (27r)4 Jo (A;2 - A(a,p2))2 ' ^ ' 

and from Eqs. (3.19), (3.21) we see that Af^ = 1. 

^^See Eq.(3.7) of ref.[ll]. The nucleon vertex function (3.11) in the noncovariant LC 
normalization is equal to the vertex function of ref.[ll] multiplied by a factor 1/ ^flpZ. 
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If we perform a partial integration in a for the term oc n — 1 in Eq.(3.19) 
and compare the resulting expression to the definition of the moments (Eq.(3.16)), 
we can directly read off the distribution function: 

(3.22) 

for < a; < 1 and zero otherwise. This expression can readily be used in Eu- 
clidean regularization schemes. Performing a Wick rotation and introducing 
the proper time regularization scheme (3.3), we finally obtain 



= - ^e--(--^) [1 + ((M. + Mf - Ml 



(3.23) 



The contribution of the diquark diagram of Fig.l to the distribution /g/Aro 
can be expressed as a convolution integral, involving the quark distribu- 
tion in a virtual diquark {fq/oo) ^-nd the diquark distribution in the nucleon 
(/n/ivo)- For simplicity we use here the on-shell approximation for the di- 
quark distribution, that is, we assume that fo/m has a sharp peak at the 
diquark virtuality go — that fq/oo can be replaced by its on-shell 

value and taken out of the integral over the virtuality. Then we use the rela- 
tion [11] j Aql fo/miy, ql) = (1 - y)' "^^cic /g^^Q was evaluated above, 
see (3.23). The contribution of the diquark diagram is then given by 

f%i^) = [ jl - y^)U/Do{z) 4/^^(1 - y) . (3.24) 

Here /^^q given by Eq.(3.23) and we anticipated the fact that fq/Do(z) 
has support for < z < 1. To evaluate fq/Do{z), we use the moments (see 
Eq.(3.10) of Ref.[ll]) 



Al?/""^ --^'9D^,j^A^\ IVz. {l,S,{k)^,So{k - q)) , (3.25) 



n-l 



dq^ J (27r)4 \q. 
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where one has to set (f = M"^ at the end. 

Before giving the exphcit expressions, let us confirm the number and 
momentum sum rules. Prom the form of go given in (3.15), it follows that 
^(9/^) _ 2^ and because we have already confirmed that the n — 1 moment 
of /q^o equal to 1, the convolution (3.24) shows that the diquark diagram 
gives a contribution equal to 2 to the number sum rule. The number sum 
rule value for the total distribution (3.17) is therefore 

^ ^xfq/No{x) =3. (3.26) 
For the momentum sum rule we obtain 

^ dxxfqlNQ{x) = ^ dxx/g^^^(x) 

(3.27) 

1 

Using the fact that fq/Do{z) is symmetric around z = - and that its n = 1 
moment is equal to 2, we see that / dz z fqiDoiz) — 1. The momentum sum 

J 

rule (3.27) then becomes 

/' dx X fq/m{x) ^ f'dxx f%{x) + f dx (1 - x) 4/^^(x) = 1 . (3.28) 

We now continue with the evaluation of (3.25) and the distribution fq/Do{x)- 
By using a Feynman parameter (a) and observing that g = we obtain 
(see Appendix D) 



^ + qM^ -a)-e + ^(1 - 2a){kl - kj) 

(A;2 + q^a{l -a)- M'^f ■ \ ■ ) 
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If we perform a partial integration in a for the term oc n — 1 in Eq.(3.29) 
and compare to the definition of the moments (Eq.(3.16), we obtain 



1 



X . o ^^7:;^^ (3-30) 



where the derivative d/dx acts on the whole function on its r.h.s. Performing 
a Wick rotation and introducing the proper time regularization (3.3), we 
finally obtain 



'-,TX[1 — X 

ml 



(3.31) 



The results (3.23) and (3.31) arc then inserted into the convolution formula 
(3.24) to give the contribution of the diquark diagram to the quark distribu- 
tion function. 

3.3 Model for nuclear matter 

The EOS of NM in the mean field approximation to the NJL model can 
be derived in a formal way [10] from the quark Lagrangian (3.2) by using 
hadronization techniques. Equivalently, but more intuitively, one can use 
the "hybrid model assumption" that the expectation value of any local quark 
operator, like the quark Hamiltonian density, can be split into the expectation 
value in the valence quark vacuum |0) = |p = 0) and an average over the 
nucleon Fermi sea consisting of correlated valence quarks. The resulting 
energy density of isospin symmetric NM at rest has the form [9, 10] 

^ = ^>'-^+V(&e(P.-|p|)%. (3.32) 
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where is given by (2.5) and the vacuum contribution (quark loop) is 

The condition d£ / dVo ~ leads to 

Fo^GG^p, (3.34) 

and we can eliminate the vector field in (3.32) in favor of the baryon den- 
sity. The resulting expression has then to be minimized with respect to the 
constituent quark mass M for fixed density. For zero density this condition 
becomes identical to the familiar gap equation of the NJL model, and for 

finite density the nonlinear M-depcndcnce of the nucleon mass A/jv, which 
is obtained from the pole of To(p) (see Eq.(3.8)), is essential to obtain satu- 
ration of the binding energy per nucleon [9] . 

4 Results and discussions 

In this Section we present our numerical results for the quark LC momentum 
distributions, the structure functions and the EMC effect. The parameters 
used in the calculation are the same as in Refs. [9, 10]. The IR cut-off Air is 
fixed at 200 MeV, while the values of m, G^^ and Auv are determined so as 
to reproduce Mq = 400 MeV by the gap equation in the vacuum, ttIt^ = 140 
MeV from the pole of the t-matrix in the pionic qq channel and = 93 MeV 
from the one quark loop diagram for pion decay. The resulting values are 
m = 16.9 MeV, = 19.60 GeV'^, and Auv = 638.5 MeV. The couphng 
constant Gg of (3.4) is determined so as to reproduce Mjvo = 940 MeV as 
the pole of the quark-diquark t-matrix at zero density (see Eq. (3.8)), which 
gives the ratio Vg — Gg/G^^ — 0.508. Finally, the couphng constant G^j is 
determined so that the curve for the NM binding energy per nucleon (Eb/A) 
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as a function of the density passes through the empirical saturation point 
(p, Eb/A) = (0.16 fm-3, 15 MeV), which gives the ratio = G^/G^ = 0.37. 
Prom the pole of the t-matrix in the vector qq channel, this value gives an uj 
meson mass in vacuum[10] of 827 MeV, which is larger than the experimental 
value (783 MeV), but is nevertheless reasonable. 

We recall from Ref.[9] that, in this simple NJL model, G^ is the only 
free parameter for the NM calculation. Therefore, although the binding 
energy curve passes through the empirical saturation point, its minimum 
is at a different point, {p,Eb/A) = (0.22 fm'^, 17.3 MeV), see Fig. 6 of 
Ref.[9]. As we discussed in Sect. 2, in the derivation of Eq. (2.39) for the 
quark distribution in NM we assumed that Mjv = cf, which holds only at 
the minimum of the binding energy curve. Therefore, although our calculated 
saturation density p = 0.22 fm~^ is larger than the empirical one, we have 
to use this density in our numerical calculations for consistency 

In Table 1 we list the effective quark, diquark and nucleon masses, the 
vector field acting on the quark ( Vq) , the nucleon Fermi energy and the mass 
per nucleon Mjv = Mm - Eb/A at p ^ 0, p ^ 0.16 fm~^, and p = 0.22 
fm~^. We note that for p — 0.22 fm~^ the sum M + Mr, is smaller than 
Mjv, which is possible because of the absence of unphysical thresholds in our 
regularization procedure. 

We first show our results for the valence up and down quark distributions 
in the free proton in Figs. 3 and 4. They are given by 



where the contributions from the quark diagram (Eq. (3.23)) and the diquark 

"'^'^However, in order to indicate the density dependence of our results, we also show in 
Appendix E the case p = 0.16 fm^-^, simply by using Eq.(2.39) for this value of the density. 



(4.1) 




(4.2) 
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p = 


p = 0.16fm-3 


p = 0.22fm-3 


M 


400 


308 


276 


Mn 


576 


413 


355 


Mm 


940 


707 


634 


Vo 





53 


75 


ep 


940 


914 


923 


Mn 


940 


925 


923 



Table 1: Effective masses M, Mjj and M^ for the quark, the diquark and the 
nucleon, the vector field for quarks (Vq), the nucleon Fermi energy (ej?), and 
the mass per nucleon Mjq = M^o — Eb/A (all in MeV) at p = 0, p = 0.16 
fm~^ and p = 0.22 fm~^. The saturation density in the present calculation 
is p = 0.22fm'~^. The values of the quark-diquark coupling constant gn (see 
Eq.(3.14)) are 18.23, 16.16, and 15.64 for p = 0, p = 0.16 fm-^ and p = 0.22 
fm~^, respectively. 

diagram (Eq. (3.24)) refer to the case p = of Table 1. 

The solid lines show our NJL model results. We see clearly that d„(x) is 
softer than Uy{x) - that is, the peak position of dy{x) is at smaller values of 
X. This difference comes from the diquark correlations [27, 12], or, in other 
words, from the difference between the "hard" distribution fg^^o 
"soft" one fq^MQ- '^'^^ total LC momentum fraction carried by the u quarks 
is 71% and that of the d quark is 29%. 

In order to make contact to the measured structure functions, we associate 
the distributions calculated in the NJL model with some low energy scale Qo, 
where a description in terms of quark degrees of freedom alone is expected 
to be vahd. We then perform the evolution according to the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations [28] up to the value of 
Q where experimental data for DIS and empirical parametrizations of parton 
distributions are avahable Here we fix Qo — 400 MeV as in our earher 
^^We use the computer code of Ref.[29] with Nf = 3, Aqcd = 250 MeV in the MS 



27 









Qo^=0.16 GeV^- 




q2-10 GeV^ ~- 




^ — ^ empirical ~ 




\ (Q^=10 GeV^) ; 


/''' / ^ 

/ / '^^ 





0.0 0.2 0.4 0.6 0.8 1.0 



Figure 3: The valence up quark LC momentum distribution in the free nu- 
cleon, multiphcd by x. The sohd hne shows the NJL model result (Qq = 0.16 
GeV^), the dashed line shows the result after the evolution up to = 10 
GeV^, and the dotted line shows the empirical parametrization of Ref.[30]. 



calculations [11, 12], and perform the evolution up to Q"^ — 10 GeV^. The 
resulting valence quark distributions are shown by the dashed lines in Figs. 
3 and 4 and they are compared with the empirical distributions of Ref. [30] 
shown by the dotted lines. We see that our calculation reproduces the main 
features of the empirical distributions. The calculated total LC momentum 
fractions at Q"^ — 10 GeV^ are 29% and 12% for the valence u quarks and 
the valence d quark, respectively, while the rest is carried by the gluons and 
sea quarks which are generated by the evolution. 

We now discuss the modifications of the quark LC momentum distribu- 
tion at finite density. For this purpose, we first replace the free quark, diquark 

renormalization and factorization scheme in the next to leading order to perform the 
evolution. 
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Figure 4: Same as Fig. 3 for the valence down quark. 



and nucleon masses by the effective ones given in the last column of Table 1, 
which gives the distribution fq/m- Then we calculate the convolution with 
the nucleon LC momentum distribution according to Eq.(2.40) to obtain 
fq/Ao- Finally we include the effect of the vector potential through the scale 
transformation (2.39) to get the total result, fg/A- The results are shown in 
Fig. 5. The dotted line shows isoscalar valence quark distribution for the 
free nucleon and the dashed line shows the result when all masses are re- 
placed by the effective ones. We see that the distribution becomes somewhat 
stiffer - i.e., it has more support for large momentum fractions. This can be 
understood as follows: For large x the spectator quark, which corresponds to 
the quark diagram of Fig. 1, gives the most important contribution. From 
the spectral representation (2.23) it follows that the typical LC momentum 
fraction carried by the spectator quark is 1 — Md/Mn - in agreement with 
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Ref. [27]. Since in our calculation the ratio Md/M^ decreases as the density 
increases (see Table 1), the fraction carried by the spectator quark increases 
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Figure 5: Sum of the valence up and down quark LC momentum distributions 
at the low energy scale (Qq = 0.16 GeV^), multiplied by the LC momentum 
fraction. The dotted line shows the distribution in the free nuclcon, the 
dashed line is obtained by replacing the free masses by the effective ones, 
the dot-dashed line shows the result of the convolution (2.40), and the solid 
line is the final result obtained from the scale transformation (2.39). The LC 
momentum fraction on the horizontal axis gcncrically refers to the proper 
variable of each distribution, i.e., x for the dotted and dashed lines, and xa 
for the dot-dashed and solid lines. 



The result of the convolution (2.40) with the nucleon momentum dis- 
tribution is shown by the dot-dashed hue in Fig. 5. Because the nucleon 

""^^Wc wish to point out that, in spite of the fact that the LC momentum distribution 
becomes stiffer, the baryon radius of the nuclcon increases in the medium [31], which shows 
that there is no direct connection between the change of the LC momentum distribution 
and the baryon radius. 
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momentum distribution fN/Ao{y) has a finite width and is symmetric around 
y — 1 (see Eq.(2.35)), the result of the convolution shows the familiar effects 
of Fermi motion. The support extends beyond x — 1 and, relative to the 
dashed line, there is a depletion for intermediate values of x. However, we 
have to note two points here. First, since at this stage the effect of the vector 
field is not yet included, the convolution (2.40) involves the nucleon momen- 
tum distribution (2.35) without the explicit effect of the vector potential, and 
has an appreciably larger width than the "full" distribution (see (2.33) and 
(2.36)). Therefore, the effect of the Fermi motion, which is seen by compar- 
ing the dashed and dot-dashed lines in Fig. 5, is more pronounced than in 
the calculations of Ref.[6], where the convolution of the quark distribution 
in a free nucleon with the full f^/A was considered. Second, Fig. 5 clearly 
shows that the EMC effect cannot be explained by taking into account only 
the effects of the scalar potential and Fermi motion: in comparison with the 
dotted line, the dot-dashed line shows an enhancement in the valence quark 
region {x ~ 0.7) instead of a quenching. 

We now come to the discussion of the effect of the mean vector field on 
the quark distribution function, as expressed by Eq.(2.39) and shown by the 
solid hue in Fig. 5. If we compare the sohd line with the dot-dashed line, 
we see that the effect of the vector field is to squeeze the quark distribution 
both from the small and large x side, leading to a reduction of the width, a 
quenching in the valence quark region and an enhancement in the region of 
smaller x. To understand the origin of this effect, we recall the key points 
which led to the scaling relation (2.39): The presence of the vector field leads 
to a shift of the quark momentum according to k'^ — kq^-V^, see Eq.(2.21) for 
H — 1, and of the nucleon chemical potential (mass per nucleon) according 
to €f — Ep + 3Vo, see Eq.(2.31). Then the ratio of the quark LC momentum 
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to the chemical potential without the effect of the vector field (x'^ 



it- 
is related to the one including the effect of the vector field {xa = — r-^) by 

which is just the rescaling of the variables expressed in Eq.(2.39). Using 
(2.31), this linear relationship between and xa can also be expressed as 

x'a ^ xa{1 + ^v) - V , (4.4) 

where v = Vo/Ep. From Eq.(4.4) we see that for = - we also have x'j^ = - 

if fq/Aoi^A) has the shape of a typical valence quark distribution with a 

peak around xa — -:z, the distribution fq/Aoi^'A — viewed as a 

3 Ep Ep 

function of xa-, will be squeezed from both sides {xa < 1/3 and xa > 1/3), 

with the amount of squeezing being larger for large xa- According to (2.39), 

ep 

this distribution is then multiplied by an overall factor — — so as to satisfy 

Ep 

the number and momentum sum rules. This mechanism, which is clearly 
seen by comparing the dot-dashed and the solid curves in Fig. 5 is largely 
responsible for the explanation of the EMC effect in a mean field description 
of NM, where the vector field couples to the quarks inside the nucleons. 

If we finally compare the dotted and solid curves of Fig. 5, we observe an 
enhancement for momentum fractions between 0.2 and 0.4, and a depletion 
for larger x, besides the enhancement for a; ^ 1 arising from the Fermi 
motion. These features are consistent with the EMC effect. 



^^The jump at small xa shown by the solid line in Fig. 5 can also be understood from 
the rescaling relation Eq.(2.39): The original distribution fq/Aoi^A) nonzero at cc^ = 
but zero for a;^ < 0, and therefore the shifted distribution fq/A{xA) is discontinuous at 
xa = Vo/ep- However, in this region one has to consider also the effect of the sea quark 
distributions, which are not included here. 
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We next show in Fig. 6 our results for the nuclear structure function per 
nucleon, calculated for isospin symmetric NM, 

5 

F2a{xa) = Y^^A fq/A^XA) , (4.5) 

in comparison to the isoscalar structure function of a free nucleon 

F2n{x) = ^ {F2p{x) + F2n{x)) = ^ X fg/mix) . (4.6) 




X 

Figure 6: The solid line shows the structure function per nucleon in sym- 
metric NM given by Eq.(4.5). It is plotted as a function of x by using the 

relation xa = — x. The dashed line shows the isoscalar free nucleon struc- 

ture function of Eq. (4.6). 



The curves of Fig. 6 were obtained by evolving the quark distributions 
shown by the solid and dotted lines in Fig. 5 up to — 10 GeV^. In order 
to plot the nuclear structure function (4.5) as a function of x, we used the 
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relation xa — — ^ x = 1.02 x, see Eqs.(2.3), (2.31) and Table 1. Although 
this "binding effect" is on the level of a few percent, it gives a non-negligible 
contribution at large x. Fig. 6 clearly shows the suppression of the structure 
function at large x and the enhancement at smaller x, which we discussed 

above for the distribution function at the low energy scale. 

From the two curves of Fig. 6, we obtain the ratio between the nuclear 
and the nucleon structure functions 



F2a{xa) 
F2n(x) ' 



RaMx) = '-^, (4.7) 



which is shown in Fig. 7. The data points in this figure are the extrapolations 
of nuclear DIS data to the NM case [32] 

We see that the calculation can reproduce the main features of the EMC 
efi^ect, namely the suppression at large x and the enhancement at smaller x. 
[The fact that in our calculation the saturation density is too large (p = 0.22 
fm"^) leads to an EMC effect that is shghtly too strong at large x. This is 
clear if we compare to the results discussed in Appendix E for p = 0.16 fm~^ 
- see Fig. 11.] 

5 Summary and conclusions 

In this paper we have discussed the modifications of the structure function 
of a nucleon bound in the nuclear medium. For this purpose, we used an 
effective chiral quark theory which can account for the quark substructure of 
a single nucleon, as well as for the saturation properties of the NM binding 
energy. Our description of NM was carried out in the mean field approxi- 
mation, which is characterized by mean scalar and vector fields coupling to 



^^The fit in Ref.[32] uses data over a wide range of Q"^, assuming that the EMC ratio 
is fairly independent of Q^. Although Fig. 7 refers to = 10 GeV^, we have confirmed 
that our results are independent of Q^, as long as Q^^4: GeV^. 
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Q^=10 GeV^ 
p=0.22 fm"^ 

Nuclear matter data: 

Sick, Day: PLB 274 (1992) 16 



0.0 0.2 0.4 0.6 0.8 1.0 

X 

Figure 7: The ratio (4.7) of the nuclear to the nucleon structure function. 
The data points are the extrapolations of nuclear DIS data to the NM case 
[32]. 

the quarks inside the nucleons. Our aim was to assess the effects which arise 
from the quark structure of the bound nucleon. 

It is well known that the effects of the mean scalar and vector fields on 
the level of nucleons, including the Fermi motion, cannot explain the EMC 
effect. We have investigated these medium modifications at the level of the 
quark structure and our findings can be summarized as follows. The mean 
scalar field can simply be incorporated into the effective masses and together 
with the Fermi motion of nucleons it cannot explain the observed reduction 
of the EMC ratio in the intermediate x region. The vector field, on the other 
hand, influences the form of the quark LC momentum distributions directly 
- on top of its indirect influence through the EOS of the system. This direct 
modification of the quark distribution by the mean vector field is expressed by 



< 

K 

0.6 
0.4 
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Eq.(2.39). As our discussion in Sect. 4 showed, it plays a key role in explaining 
the EMC effect, in the framework of a mean field description of NM in terms 
of quarks. We note that the relation (2.39) includes the combined effect of the 
mean vector field on the quark momentum distribution in the nucleon and 
on the nucleon momentum distribution in NM. While the effect of the mean 
vector field on the nucleon momentum distribution alone tends to cancel the 
effect of the scalar field [6], its combined effect is appreciable and works in 
a direction which is consistent with the EMC data. It is therefore crucial 
to include the quark substructure of the nucleon for the explanation of the 
EMC effect. 

Our findings are reminiscent of the discussions in Ref. [9] (and, of course 
Refs. [14]), which emphasized the importance of the fact that the mean scalar 
field couples to the quarks in the nucleon, instead of to the nucleon as a 
whole, for the stability of the NM ground state in chiral effective quark 
theories. Concerning the mean vector field, on the other hand, the quark 
substructure of the nucleon is not essential as long as we consider only the 
bulk properties of the system. It is essential, however, for the structure 
functions and the EMC effect, as we have shown in this paper. The further 
study of the influence of the quark substructure on medium modifications of 
nucleon properties along these lines is an important aspect of what has come 
to be called quark nuclear physics. 
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Appendices 



A Verification of Eq.(2.18) from Feynman di- 
agrams 

In this Appendix we show that the result of an exphcit evaluation of the 
Feynman diagrams shown in Fig.l satisfies the general relation (2.18). 

First we consider the quark diagram of Fig.l, which gives the following 
contribution to the distribution function (see Eq.(3.7) of Ref.[ll]): 



ft/Ni^) = r^(p) [^^n{x.p)\ Vn{p) (A.l) 



with 



niv(x,p) = - / ^5 [x - - j S{k)T{p - k) . (A.2) 

Here we make use of (see (3.7) and (3.5)) 

S{k) ^ S^{kc^) ■ T{p-k)^To{{p-k)^)^T^{p^-kq), (A.3) 

where the kinetic momenta are defined as in (2.21) for H=Q, D and N, 
corresponding to the quark numbers = 1, 2, 3. If we make a shift k — > 
k-\-V \n (A.2), the 5 function becomes 

V P- J Pn- V Pn-/ 

p- V_ 

where x' = x . By noting that the on-shell value of p% is Ep, see 

Pn- Pn- 

Eq.(2.5), and that the nucleon spinor rjv(p) of Eq.(3.11) differs from the one 
for a free nucleon (quark-diquark bound state) only by the replacement of the 
free nucleon mass by the effective one, we immediately obtain the relation 

' Pn- ' Pn- Pn- 
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where the distribution /q/7vo(^') differs from the expression for a free nucleon 
only by the replacement of the free masses by the effective ones. 

In the on-shell approximation for the diquark, to which our actual calcu- 
lations of Subsect.3.2 refer, the contribution of the diquark diagram of Fig.l 
to the quark distribution function has the form [11] 

^ j j '^(^ ~ yz)fq/D{z) fD/N{y) , (A.6) 

where 

fn/Niy) = 4/^^(1 -y) = ^fn/m (y' -Z^y-^)- (A-7) 
' Pn- \ Pn- Pn-/ 

To derive the second equality in (A.7), we used the result (A. 5), and fo/miy') - 
/i/7v(l ~ y')- '^^^ quark distribution in an on-shell diquark is given by (cf. 
Eq. (3.19) of ref.[ll]) 

at the on-shell value — ^d- Here the scalar bubble graph is given by 
(3.6), and 

n.(^, q)^6tj ^5 (^z - tr,, [75 S{kh, S{-{q - k))] . (A.9) 

Here we use again the relations S{k) = 5'o(A;q), S{—{q — k)) = So{—{q — 
^)q) — So{kQ — q'd). Then we perform a shift k —>■ k + V , and rewrite the 
delta function as 

M, (A.10) 

q- V- 

where z = z . Noting also that Iis{q) — 1150(^0), we obtain 

fq/D{z) = — fq/DO (z = —Z - — ) . (A.ll) 
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The distributions (A. 7) and (A. 11) must then be inserted into the convolution 
(A. 6). Here we must note that, because of the presence of the vector field, 
the distribution (A. 11) actually depends not only on but only on or on 
y (since y — q^/p^). Convolution integrals of this type are treated generally 
in Appendix C, and here we have a special case of (C.l) ior h — D and 
H — N. Therefore we can immediately read off the result from (C.ll): 

where the distribution fg^^oi^') differs from the free nucleon expression only 
by the replacement of the free masses by the effective ones. The sum of (A. 5) 
and (A. 12) then reproduces the more general relation (2.18) of the main text. 

B Spectral representations of LC momentum 
distributions 

In this Appendix we show that the support of fq/niz) given by Eq.(2.22) is 
consistent with the spectral representation (2.23). Following the discussion 
of Ref.[21], the contributions to the spectral sum on the r.h.s. of (2.23) can 
be represented by the four diagrams of Fig. 8 
In all four diagrams, we have from (2.23) 

p'l^ = p.{l - z) . (B.l) 

In the following we discuss the support from these diagrams separately: 

1. Diagram 1: The intermediate state \n > has quark number {H — 1), 
where H is the quark number of the initial hadron H. Because the 



""^^In general, the arrows in these diagrams denote the flow of quark (or baryon) number: 
For example, the left vertex in diagram 1 means to remove a quark with momentum 
P- X, or to add an antiquark with momentum —p- x, and vice versa for the right vertex. 
However, since the calculations of the main part are restricted to the valence quark picture, 
we need not consider the aniquarks here. 
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1 



2 



P_V\ \/P-2 



P_V\ \/P-2 




p_z/\ \/p_z 




P 



P-V\ 




\/p_z 



Figure 8: Diagrammatic representation of the r.h.s. of the spectral represen- 
tation (2.23). The bold line denotes the external hadron H, the thin line a 
quark, and the wavy line the intermediate state. 



minus component of the kinetic momentum must be positive definite, 

we obtain 

-(H- = - - (// - > , 
which gives the support 

z<l- (-^ ~ 1)^- (diagram 1) . (B.2) 
P- 

2. Diagrams 2 and 3: The intermediate state is \n >— \n' > <Si\H,p >, 
where \n' > is indicated by the wavy line in the diagrams, and has 
quark number —1 and minus momentum component given by 

(n') (n) /T^ o\ 

p_ — — P- — —p- Z , 

where we used (B.l) in the second step. Since the minus component of 
the kinetic momentum of the state \n' > must be positive definite, we 
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obtain 



which gives the support 

z < — (diagrams 2 and 3) . (B.4) 

P- 

3. Diagram 4: The intermediate state is \n >— \n' > <Si\H,p > <Si\H,p >, 
where \n' > is indicated by the wavy hne in the diagram, and has quark 
number —{H + 1) and minus momentum component given by 

pin') ^ pin) _ ^ _p_ ^ 

where we used (B.l) in the second step. Since the minus component of 
the kinetic momentum of the state \n' > must be positive definite, we 
obtain 

pin') + (^H + 1)V^^ -p_ {l + z) + {H + 1)V. > , 
which gives the support 

z<-l+ + (diagram 4) (B.6) 

P- 

We therefore see '^^ that in the region 

^<,<:_(^L^ (B.7) 

only diagram 1 contributes. Since diagram 1 corresponds to the "quark 
distribution" in the parton model sense, we obtain the support as given by 
Eq. (2.19). 

^°Note that the r.h.s. of (B.6) is smaller than because the minus component of the 

P- 

kinetic momentum of the external hadron H is positive: p- — HV- > 0. 
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C Convolution integrals 



Here we consider convolution integrals of the form 

fq/H{x) ^ j dy J dz S{x - yz) fq/h{z) fh/H{y) , (C.l) 

where the distributions refer to hadrons H and h with quark numbers H and 
h, respectively, and the dependence on the mean vector field is expressed as 

p / p hV_\ 

fh/H{y) = fh/m \y' = -^y - 1 (C.2) 

fq/h{z) = — fq/h0\z' = —Z-—\ . (C.3) 

The minus momentum components of the hadronic states H and h are P_ 
and p_, respectively, and their kinetic momenta are defined as in (2.7). For 
the case of the convolution (2.8), which is needed in the main text, we have 
to identify H=A, h=N, and P = ep/V^, Pr^ = Ep/y/2. 

To evaluate (C.l), we make the transformation of variables 

{y,z) {p-,k-), where 

and obtain from the Jacobian of this transformation 

dy dz = — ^dp_ dA;_ . (C.5) 

P- p- 

In terms of these variables we have x — k^/P^, and the quantities y' and z' 
in (C.2) and (C.3) are expressed as 

y' = ^ z' = ^. (C.6) 

Ph- Ph- 



Then Eq.(C.l) becomes 

fq/H(x) = Jdk^l dp_ s(x-'^^ p^J^'^' ( t;:) ^^"-^ 
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j dkq- j dph- 



S X- 



ko- + V. \ 1 



where we performed a shift {k-,p-) {kQ-,ph-) in the second step. We 
then introduce the variables 

-Ph- Ph- 

p^_^Pu-y', kQ^^Pn-y'z' (C.8) 

V Ph- + HV_J Ph- V^h- Ph-J 
^- 5{x'-z'y') , (C.9) 



where in the last step we defined 



P V 

x'^-^x--^. (C.IO) 

-TH- -TH- 



Then we obtain finally 

P 



Uih{x) = ^Jdy'Jdz'S{x'-y'z')f,/Ho{y')f,/ho{z') 
^ ( P V \ 

-U/HOW = -—X--—] . (C.ll) 
H — \ H — H — / 



P, 



D Moments and distribution functions in the 
proper time regular izat ion scheme 

In this Appendix we derive the quark distribution functions in the absence 
of the vector field given in Subsect.3.2. 

The contribution of the quark diagram of Fig.l to the moments is given 
by (3.18), where only the pole term of Tq (see Eq.(3.14)) contributes. We 
combine the two denominators by a Feynman parameter a, and shift the 
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loop momentum according to k ^ k + pa to obtain 

where A(q;,p^) is defined in Eq.(3.20). In general, invariant integration over 
a function of the form f{k^)k"'k^ . . . gives a sum of products of contractions 
g"'^ . . .. If we note that g = 0, it is easy to see that, in the expansion 
of {k- +p_Q;)"~^, all terms oc (k^)"^ give no contributions if m > 2. We 
can therefore replace +p_Q;)"~^ —>■ (p_Q;)"~^ + (n — (p-a)""^. If 
we further use — /i;+7_ + A;-7+ — fcj. • 7j^, as well as the form (3.11) of the 
nucleon spinor, we obtain 



P- 

n-l 



— / TTT^ / da a 
dp+ J 27r 4 io 



a^+M + ^{kl-kl) 7 
X — 



dp+J (27r) 

un{p) ■ 



{k-' - A{a,p'')f 

By using the Dirac equation, we can replace 7_ p_lMj^ in the last term 
of the above expression. Because of p^ = the integrand can then be 
considered as a function and by noting that 

we arrive at the form (3.19) given in the main text. If we perform a partial 
integration in a and note that the surface term is independent of we obtain 
the expression (3.22) for the distribution function. 

We perform a Wick rotation in (3.22), and note that in any 0(4) invariant 
Euchdean regularization scheme we can replace k^ + k^ ^ ^^/2, where now 
our notation refers to Euchdean metric (A;^ — k^ + k"^). By introducing the 
proper time regularization scheme (3.3), we obtain 
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It is then easy to perform the loop integral and the required derivatives to 
arrive at the result (3.23) of the main text. 

Next we consider the quark LC momentum distribution in the diquark. 
The moments are given by Eq.(3.25). We calculate the Dirac trace, introduce 
a Feynman parameter a, and shift the loop momentum according to A; — > 
k + qa to obtain 



where one has to set — M"^ at the end. Following the same argu- 
ment as given above for the quark diagram, there are no contributions from 
(m > 2) in the expansion of (/c_ + g_Q;)"~^ because of g = 0. By 
using k ■ q = k^q^ + k^q^ — k± ■ q^, we obtain the expression (3.29) of the 
main text. Performing a partial integration in a and noting that the surface 
term is independent of q^ leads to the distribution function (3.30). We then 
perform a Wick rotation in (3.30), and note again that in any 0(4) invariant 
Euclidean regularization scheme we can replace + ^ ^^/2, where now 
k^ — kg + k^. By introducing the proper time regularization scheme (3.3), 
we obtain 



where the derivative d/dx acts on the whole function on its r.h.s. It is then 
easy to perform the loop integral and the required derivatives to arrive at 
the result (3.31) of the main text. 




X 



+ q'^ajl -a)-e + k- q{l - 2a) 
(A;2 + q'^a{l -a)- M^f 
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E Results for the case p =0.16 fm 



The numerical calculations discussed in Sect. 4 were performed at the satu- 
ration density in our model (p = 0.22 fm""^), which is too high compared to 
the empirical value. In order to indicate the dependence on the density, we 
will show the results for p = 0.16 fm~^ in this Appendix, using the effective 
masses and energies listed in Table 1. One should keep in mind, however, 
that in the derivations of the main text, for example of Eqs.(2.34) and (2.39), 
the relation = ep was assumed. As this is satisfied only at the saturation 
density (see Table 1), the results shown in this Appendix should serve only 
to indicate qualitatively the density dependence of the results. 

Fig. 9 shows how the results of Fig. 5 change if we decrease the density 
to p = 0.16 fm-3. 




0.00 0.25 0.50 0.75 1.00 

LC momentum fraction of quark 
Figure 9: Same as Fig.5 for the density p — 0.16 fm~^. 
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We see that the various medium modifications of the quark distribution 
function show a very similar behavior, but their magnitudes are smaller be- 
cause of the smaller density. The structure function per nucleon in nuclear 
matter (see Eq.(4.5)) is compared to the free nucleon structure function in 
Fig. 10. In order to plot the function (4.5) against x, we used the relation 
(2.3) with Mat equal to the mass per nucleon at the density p — 0.16 fm~^ 
given in Table 1. 




The quahtative behavior is the same as in Fig. 6, but the difference 
between the dashed and sohd fines is smaller. Finally, tfie ratio (4.7) is 
shown in Fig. 11. Compared to Fig. 7, tfie magnitude of tfie EMC effect is 
reduced. 
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